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Jet quenching effect has been investigated in the direct photon production, based on a realistic 
data-constrained (3+1) dimensional hydrodynamic description of the expanding hot and dense mat- 
ter, a reasonable treatment of the propagation of partons and their energy loss in the fluid, and a 
systematic study of the main sources of direct photons. Our resultant pt spectra agree with recent 
PHENIX data in a broad pt range. Parton energy loss in the plasma eventually effect signiflcantly 
direct photon production from fragmentation and jet photon conversion, similar to hadron suppres- 
sion in central heavy ion collisions. But this only causes about 40% decrease in the total production 
of direct photons, due to the mixture with other direct photon sources. 
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I. INTRODUCTION 

The formation and observation of a quark-gluon 
plasma (QGP) in heavy ion collisions are important goals 
of modern nuclear physics P, Q] . Suppression of high pt 
hadron yields 01 is one of the most important features 
observed at the Relativistic Heavy Ion Collider (RHIC). 
Theoretically this is attributed to the interaction between 
jets (hard partons) and the bulk matter [1, H, This 
offers us an excellent opportunity to study the interaction 
of partons inside the system and, consequently, proper- 
ties of the matter under extreme conditions. 

However the observed centrality-dependence of 
hadronic jet quenching doesn't show up in direct photon 
production^. In this paper, we investigate the role 
of jet quenching on direct photon production. For 
this purpose, we need a careful consideration of three 
components: 1) the space-time evolution of the hot dense 
matter or the bulk matter; 2) hard parton production 
and propagation in the bulk and 3) all sources of direct 
photons. 

The space-time evolution of the hot dense matter is 
achieved by using three dimensional (3D) hydrodynamic 
simulations of bulk matter [lH which have already 
been tested against a vast body of low pt hadron data at 
RHIC. 

The interaction between hard parton and the bulk so 
that energy loss of the hard parton is formulated via the 
BDMPS framework Q and tested on pion suppression at 
various centralities. 

Four significant sources are considered in this paper: 
leading order contribution in Primordial NN scattering, 
thermal contribution, contribution from jet-photon con- 
version and from jet fragmentation. There are pos- 
sible sources ignored at present study: The medium- 



induced photon radiation is expected to contribute at low 
transverse momentum region where thermal contribution 
dominants. The interaction involving non-equilibrated 
soft partons may also produce direct photons, but the 
preequilibrium time interval is much shorter than the life 
time of the equilibrated matter 20 fm/c). 

Therefore energy loss of hard parton in the bulk will 
directly effect two direct photon sources, jet-photon con- 
version and from jet fragmentation. We will see the role 
of jet quenching on direct photon production clearly via 
investigating those production in Au-Au collisions at dif- 
ferent centralities where various sizes of hot dense matter 
are formed. Certainly different sizes of bulks will also 
contribute thermal photon production differently which 
will not be discussed in detail in this paper. 

The paper is organized as follows: In Sec. m we first 
give a brief review on the space-time evolution of the hot 
matter created in Au+Au collisions at different centrali- 
ties based on a (3+l)-dimensional ideal hydrodynamical 
calculation. In Sec. Illlj we discuss parton energy loss 
in the QGP. We investigate neutral pion production in 
the high pt region in order to fix the parameters of the 
energy loss scheme. We discuss sequently the contribu- 
tions from various sources to direct photon pt spectra in 
Sec.lIVI We show our results and compare them with re- 
cent experimental data in Sec.[Vl Section [Vl] is devoted 
to conclusion of the present study. 



II. SPACE-TIME EVOLUTION OF THE HOT 
AND DENSE MATTER 

In our calculation, fully three-dimensional (3D) ideal 
hydrodynamics [lfl,[lH is employed to describe the space- 
time evolution of the hot and dense matter created in 
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Au+Au collisions at RHIC energy at various centralities. 
We solve the equations of energy-momentum conserva- 
tion 



d^Tt"" = 



(1) 



in full 3D space {T^x,y,rf) under the assumption that 
the local thermal equilibrium is reached (maintained) at 
(after) an initial time tq =0.6 fm/c. Here r and 77 are the 
proper time and the space-time rapidity, respectively, x 
and y are transverse coordinates. In the transverse plane, 
the centers of two colliding nuclei are located at (x, y) = 
(&/2, 0) and (—6/2,0) before the collision at an impact 
parameter b. Ideal hydrodynamics is characterized by 
the energy-momentum tensor. 



T^"' = (e + P)u^'u" - Pg^"', 



(2) 



where e, P, and are energy density, pressure, and lo- 
cal four velocity, respectively. We neglect the finite net- 
baryon density which is small near the mid-rapidity at 
RHIC. For the high temperature (T > = 170 MeV) 
QGP phase we use the equation of state (EOS) of mass- 
less non- interacting parton gas (u, d, s quarks and glu- 
ons) with a bag pressure B: 



p= i(e-4B). 



(3) 



The bag constant is tuned to be Bi = 247.19 MeV to 
match pressure of the QGP phase to that of a hadron 
resonance gas at critical temperature Tc = 170 MeV. 
A hadron resonance gas model at T < Tc includes all 
hadrons up to the mass of the A(1232) resonance. Our 
hadron resonance gas EOS implements chemical freeze- 
out at Trh = Tc = 170 MeV, as observed in collisions at 

RHIC [Hi. 

We assume that, at tq = 0.6 fm/c, the initial entropy 
distributions is proportional to a linear combination of 
the number density of participants (85%) and binary 
collisions (15%) Centrality dependence of charged 

particle multiplicity observed by PHOBOS has been 
well reproduced by full 3D hydrodynamics simulations 
with the above setups [l3|. In the following calculations, 
hydrodynamic outputs at representative impact parame- 
ters b = 3.2, 5.5, 7.2, 8.5, 9.7, and 10.8 fm are chosen for 
0-10%, 10-20%, • • • , 50-60% centrality, respectively. 

For convenience of the following calculations, we intro- 
duce /qgp(''') X, y, rf) as the fraction of the QGP phase in 
a fluid element. It is obvious that /qgp = 1 (0) in the 
QGP (hadronic) phase. In the mixed phase, the fraction 
of the QGP is calculated via 

/qgp eQGP + (1 - ./QGp)ehad e(T, x, y, 77) 

with bqgp and Chad being the energy densities of the 
QGP phase and the hadron phase at T = Tc, respectively. 



III. PARTON ENERGY LOSS IN A PLASMA 

Hard parton production from primordial nucleus- 
nucleus scattering is calculated as [l5l | 

= KTAB{b) J dXadXbGa/A{Xa,M^) 



abed 



xGusixb, M^)-^iab ^ cd)6{s + i + u) (4) 

TT dt 

where TAsib) is the nuclear overlapping function at an 
impact parameter b for each centrality, Ga/A{xa, M^) and 
Gb/B{xb, M'^) are parton distribution functions in nuclei 
A and B. We take MRST 2001 LO parton distributions 
in proton [lB|. if = 2 is chosen to take into account 
higher order contributions. These parameters are chosen 
as to reproduce high pt pion data in pp collisions at RHIC, 
which will be discussed later. Nuclear shadowing effect 
and EMC effect are taken into account thr oug h EKS98 
scale dependent nuclear ratios R^^^{x,A) [l3|. Isospin 
of a nucleus with mass A, neutron number N, and proton 
number Z is corrected as follows: 



Ga/A{x) 



-jGa/N{x) + -^Ga/p[x) 



R 



EKS 



{x,A). (5) 



The isospin mixture and nuclear shadowing eventually 
cause a decrease of nuclear modiflcation at high pt region. 

We assume that all jets are produced atT = l/(5«0 
with the phase space distribution 

/o(p, r) oc '^Ta [x - ^, y^ Tb + |, ^{z) (6) 

where f = {x,y,z) is the coordinate of a jet, b is the 
impact parameter, and Ta and Tb are thickness functions 
of nuclei A and B. 

When hard partons propagate in a hot dense matter, 
they may lose energy in the medium via radiation or 
elastic collisions with thermal partons0, 0, @, 0I- This 
is a quite hot issue and several different schemes have 
been employed for jet energy loss in QGP phase while 
the study of hard parton energy loss in hadronic phase 
is not so popular. In this paper the BDMPS schemes is 
employed which considers the energy loss of hard parton 
in QGP phase via radiation. A free parameter D is intro- 
duced by us to take into account of energy loss via other 
modes. 

The total energy loss along the trajectory of a hard 
parton is calculated as 

AE{t,po,ro) = D / dTe{i,T,x{T))e{fQGp{T,x{T))), 

J To 

(7) 

where x{t) is the trajectory of parton, /qgp(t, a;(r)) is 
the fraction of the QGP phase at a position {t,x{t)), and 
6* is a step function, which gives 6'(/qgp) equal unity in 
the QGP and the mixed phases and zero in the hadron 
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phase. In the mixed phase, the thermal parton den- 
sity will be diluted with the factor /qgp, the fraction 
of the QGP phase. Here D is an adjustable parameter, 
e(i, T, x{t)) is the energy loss per unit distance for a par- 
ton z at a position {t,x{t)), given as 



e(i, r, x{t)) = as\J [i^E* jXi. 

Note that the above quantities, i.e., temperature T, fluid 
velocity w^i, parton densities Pi and, in turn, mean free 
path Ai, depend on the location of the parton x{t) and 
can be obtained from full 3D hydrodynamics simulations 
discussed in the previous section. 

We first discuss pion production in proton-proton col- 
lisions. We calculate neutral 7r-meson production assum- 
ing pQCD factorization, Eq. (|4]), 



dyd'^pt 



J2 / 



dNPP- 



1 



dyd'^p^ 



(8) 



where D'^o/ci^cQ'^) is 7r° fragmentation functions pa- 
rameterized by Kniehl et al. (l8l |. In the high pt region 
where the pQCD is expected to work, we reasonably re- 
produce the experimental data with the above setup with 
K = 2 and M = Q = pt. We use the pt spectrum as a 
reference spectrum in the following calculations. 

The effect of parton energy loss is taken into ac- 
count through the medium modified fragmentation func- 
tion [Hi D^o/^{zc,Q'^,AEc). With a common value of 
the parameter D = 1.5, we can reasonably reproduce the 
PHENIX data [§| on 7r° jet quenching in the high pt re- 
gion at all centralities simultaneously. In the following 
photon calculations, we always use the BDMPS energy 
loss formula (d with D = 1.5. 



IV. THE DIFFERENT SOURCES OF DIRECT 
PHOTON PRODUCTION 



Possible sources of direct photons are considered in the 
following. 

Leading order contribution in Primordial NN scatter- 
ing. The direct photon production via Compton scatter- 
ing and quark-antiquark annihilation can be calculated in 
perturbation theory using the conventional parton distri- 
bution functions and the factorization hypothesis. Higher 
order contributions are considered as a part of the so- 
called jet fragmentation contribution and jet photon con- 
version which will be discussed separately. Similar to 
Eq. (jU, the leading order contribution to direct photon 
production in nucleus-nucleus collisions reads 

— ^-^j— = Tab (&) ^ J dxadxbGa/A{xa,M'^) 
ab 

xGb/B(xfc, M^)-^{ab ^ 7 + X)S{s + t + u) (9) 

TT dt 

where the elementary processes ab j + X are Compton 
scattering qg jq and annihilation qq gj. 



Thermal production. In high energy nuclear collisions, 
the density of secondary partons is so high that the 
quarks and gluons rescatter and eventually thermalize 
to form a bubble of hot QGP. Thermal photons can be 
produced via interaction between bulk particles during 
the whole history of the evolution of hot matter. The 
interaction of bulk particles and the bulk properties can 
be refiected by the emission rate of direct photons in 
the QGP phase[2l| and in the hadronic phasefl^l, noted 
as r = Ed^R/d^p, where R is the number of photons 
emitted from a medium per unit space-time volume with 
temperature T . Total yields of thermal photons can be 
obtained by summing the emission rate over the space- 
time volume as 



dN' 



thermal 



dyd^Pt 



J d^xr{p'^u^,T) (10) 



with d'^x = Tdrdxdydri, the photon's four momentum 
in the laboratory frame and a local fiuid velocity. 

Jet fragmentation. At any stages of the evolution of 
a jet (final state parton emission), there is a possibility 
of emitting photons. The existence of a QGP affects the 
results because energetic partons lose their energy prior 
to fragmentation. In this work, we assume fragmentation 
of partons only outside the plasma, which is similar to 
high Pt hadron production from jet fragmentation. 

In pp collisions, there is no formation of QGP and no 
energy loss of hard parton in QGP. Higher order contri- 
butions to direct photon production can be considered as 
jet fragmentation and be calculated as 



dyd^Pt 



E 



dZr 



dNPP- 



1 



dyd^pl 



D'^j^z^^Q^), (11) 



with photon fragmentation functions D°y^(2;, Q^) being 
the probability for obtaining a photon from a parton c 
which carries a fraction z of the parton's momentum. So 
Pt — Pt/zc is the transverse momentum carried by the 
parton c before fragmentation and d^p/E = z^d^p'^/E'^. 
The effective fragmentation functions for obtaining pho- 
tons from partons can be calculated perturbatively. We 
use the parameterized solutions by Owens [l5| . 

In case of heavy ion collisions where a QGP is formed, 
energy loss of hard partons in a plasma should be taken 
into account. This can be done via modified fragmenta- 
tion functions [191 



D^/,(Ze,0',ASe) 



= l-e 



^°/c(^c,g'), 



(12) 



with z'^ = pt/ipl - Ai^c) and z'g = (L/Ac)pt/A£'c being 
the rescaled momentum fractions carried by the parton 
c and the emitted gluons before fragmentation. Ac is 
mean free path of the parton c in the plasma, L is the 
path length of each parton traversing the plasma. Thus, 
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in heavy ion collisions, contributions from fragmentation 
become 



dN 



AB 



dycP-pt 



E 



dz, 



dN 



AB- 



1 



(13) 

Jet photon conversion with jet energy loss. When hard 
partons propagate in a plasma, they also collide with 
thermal partons and produce direct photons via Comp- 
ton process and the quark-antiquark annihilation. We 
call this process jet-photon conversion, since it is a con- 
version of a jet into a photon with almost the same mo- 
mentum as the one of originating jet parton. Contribu- 
tion from the jet-photon conversion is calculated by inte- 
gration of conversion rate over the space-time evolution 
of the hot and dense matter in the QGP phase 



dN^P'= 
dyd^pt 



pp<^{E*,T)fQGpi^,y,V.r)d^ 



(14) 



The photon production rate by annihilation and Comp- 
ton scattering of hard partons in the medium [13, [23 | 



47r2 



q 



e-gfq{p,x)T 



■ 4^2^ 
In^ C 



(15) 

where E* is the photon energy in the local rest frame, 
C = 2.323, m^j^ = and the strong coupling 

as = (7^/47r being temperature dependent as in Eq. (??). 
a is the electromagnetic couplings, Cq and fq{p,x) are 
the electric charge and the phase-space density of a hard 
parton of flavor q. The phase space distribution of hard 
partons at t is obtained by considering parton energy 
loss as 

f{p,x) = f{p,f,T) 

= / d^PafoiPo, r- vt)S{po - p- vAE) 



where /o {p, r) is the phase space distribution at t = 
described in Eqs. The (S-function expression reflects 
the energy loss of a parton moving along a straight line 
trajectory with v = p/E = po/Eq. AE is the energy 
loss from tq to t and calculated similar to Eq. ^ but 
replacing the upper limit of integral oo with r. 



V. RESULTS AND DISCUSSION 

In Fig. [U the calculated pt spectra of direct photons 
in Au+Au collisions at y^J/vjv = 200 GeV at centrality 
0-20% and 20-40% are compared to PHENIX data 
[2^ . The PHENIX data are reproduced within our multi- 
component model remarkably well. 

The nuclear modiflcation factor Raa is obtained by 
dividing a pt spectrum in nucleus-nucleus collisions by 
the -/Vcoii-scaled pt spectrum in pp collisions. The calcu- 
lation of direct photons in pp collisions includes the lead- 
ing order contribution plus fragmentation contribution. 



In high Pt regions, our result agrees with the PHENIX 
data[2a|j|20| reasonably well: So we use it to calculate nu- 
clear modification. It also provides a baseline calculation 
with the LO contribution and fragmentation contribution 
in Au+Au collisions. 

Figure [2] shows how the nuclear modification factor for 
direct photons, Raa, depends on centrality and on energy 
loss. Data for 0-10% centrality are taken from Refs. [23| 
and [2^ . Figure[2ja) shows centrality dependence of Raa 
compared to the PHENIX data. The three curves are 
respectively 0-10% (dotted line), 20-30% (soHd line), and 
40-50% (dash-dotted line). 

Raa has a weak centrality dependence at high pt re- 
gion. This result is consistent with the observed phe- 
nomenon [i^l that the pt-integrated (for pt > 6 GeV/c) 
Raa of direct photons is almost independent of collision 
centrality. Does this imply a very weak effect from jet 
quenching? Figure [2ljb) answers this question (here for 
the most central collisions) : Comparing calculations with 
(dotted line) and without energy loss (solid line), one 
finds a difference of up to 40%. So the effect of parton en- 
ergy loss is quite visible in the pt range between 4 GeV/c 
and more than 20 GeV / c. If we would do the Raa calcu- 
lations without energy loss, the difference between cen- 
tral and semiperipheral collisions would be about 20%, 
wheras the complete calculation gives the same result 
for all centralities, within 5%. We have to admit that 
we talk about small effects, requiring experimental data 
with relative errors of less than 5to observe the effects. 

Parton energy loss in the plasma suppresses the frag- 
mentation contributions and jet-photon conversion. So 
we study the ratios of the contribution with energy loss 
to the one without energy loss, as shown in Fig. [3] ((a) 
for fragmentation and (b) for jet-photon conversion). En- 
ergy loss in the plasma depends on the path length of the 
hard parton inside the plasma, which turns out to depend 
on the collision centrality. We do see a similar centrality 
dependence of the suppression for 7r° (jet quenching ef- 
fect) in fragmentation contributions and jet-photon con- 
version. 



VI. CONCLUSION 

We calculated the centrality-dependence of pt spectra 
for direct photons in Au+Au collisions at the RHIC en- 
ergy, based on a realistic data-constrained (3+l)-D hy- 
drodynamic description of the expanding hot and dense 
matter, a reasonable treatment of propagation of par- 
tons and their energy loss, and a systematic considera- 
tion of main sources of direct photons. In this study, 
four main sources are considered, namely, leading order 
(LO) contribution from primordial elementary scatter- 
ings, thermal radiation from the fiuids, fragmentation 
from hard partons, and jet photon conversion (JPC). 
Similar work [2§| has been done before the appearance of 
the most recent data [2^. Our results agree nicely with 
the recent low pt data. 



5 



Au+Au^Y Vsfj^=200GeV 0-20% 



. PHENIX08 
o PHENIX07 




20 

p, (GeV/c) 




20 

P( (GeV/c) 



Figure 1: (Color Online) Direct photon production in Au+Au collisions at centrality 0-20% and 20-40%. PHENIX data are 
shown as open circles [2J| and filled circles (26l |. 
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Figure 2: (Color Online) The nuclear modification factor of direct photons in Au+Au collisions Raa- Data for 0-10% centrality 
are from Ref. ^ and Ref. [H]. (a): Raa at centrality 0-10%) (dotted line), 20-30% (solid line), and 40-50% (dash-dotted line) 
respectively, (b): Raa at 0-10% centrality with energy loss (dotted line) and without energy loss (solid line). 



The role of jet quenching in the high pt region of direct 
photons production has been checked via fragmentation 
photons and jet photon conversion sources. For these 
two sources, the suppression of the photon rate due to 
parton energy loss is significant in central Au+Au colli- 
sions, and becomes less important towards peripheral col- 
lisions, similar to the suppression for meson production. 
Since experimentally one may separate isolated photons 
(LO+JPC) and associate photons (fragmentation pho- 
tons), our prediction may be tested in the future. 

Considering the total yields of direct photons, the con- 
tribution from fragmentation and conversion are small, 
contributing between 5% and 10%. However, parton en- 
ergy loss plays nevertheless an important role: Without 
it, these second order effects would contribute signifi- 
cantly. Without jet quenching, the nuclear modification 



factors Raa would depend visibly on the centrality of 
the collisions. A strong energy loss is actually necessary 
to get the centrality scaling of Raa in our calculation - 
a centrality scaling which has observed by the PHENIX 
collaboration. In this sense, properties of the bulk mat- 
ter affect the photon yields at intermediate values of pt , 
via the parton energy loss. 
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